Non-perturbative Effects Near the Upper Critical Field in Type— II Superconductors 
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The properties of a conventional type-II superconductor near the upper critical field are usually 
calculated using a perturbative expansion in the strength of the order parameter. Here we show 
that perturbation theory breaks down near the upper critical field and that this breakdown leads to 
unusual behavior near the phase transition. 
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Standard microscopic calculations of the properties of 
type-II superconductors near the upper critical field rely 
on two assumptions: 1) a perturbative expansion in the 
strength of the superconducting order parameter and 2) 
weak magnetic field. These calculations yield an upper 
critical field H C 2 (T) close to what is observed experimen- 
tally jl]^] . More complete solutions to the Gor'kov equa- 
tions, which do not rely on the weak magnetic field as- 
sumption, yield low temperature divergences in H C 2(T) 
|j|-|5| which have not been observed. In this paper we 
show that the origin of these divergences is the break- 
down of perturbation theory, and that the full, non- 
perturbative solution to the Gor'kov equations describes 
tails of residual superconductivity above H C 2 . 

The microscopic description of type II superconductors 
is based on Gor'kov's extension of the BCS theory 



ft = J dr ct a H (r)c ra - -V Q J dr^ta' 1 



(1) 



which describes a short-range attractive ineraction of 
strength Vq > for particles near the Fermi energy, 
Ho is the bare Hamiltonian Hq = (l/2m) 
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(iV - eA/cf - E F , c\ a are the electron creation oper- 
ators, and fl is the system volume. The variational ap- 
proach to this problem is to solve for the eigenstates of 
the two-body Hamiltonian 



where 
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'~ri L ri J • Minimizing (Ti) leads to the self- 
consistency condition 
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The standard result for the upper critical field H C 2 follows 
from evaulating Eq. (||) perturbatively in <&(!■). Both 
sides of this equation are linear to lowest order so set- 
ting the resulting coefficients equal gives the condition 
which determines H C 2- Although this is usually done us- 
ing Green's functions, here we consider the wavefunctions 
directly. 

First we separate the magnitude of the order parameter 



$(r) = </>/(r) 
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where <f> is real and positive and /(r) is normalized so 
that ^ J |/(r)| 2 dr = 1. Next, we write the eigenstates of 
the bare Hamiltonian as 



Ho ip a (r) = £ a tpa(r) 



(•5) 



and define new quasiparticle operators by rotating to this 
basis: d* aa = Jdr ^i ce (r)c4 <T . We can then define a pairing 
matrix 
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so that the variational two-body Hamiltonian is 



H' = ^ 
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where ^ is the vector of quasiparticle operators ^ = 

[ • • • d * do, • • • ] and e is the diagonal matrix of 

bare eigenvalues |e a | < E c . 

Eq. (Q) is the general form of the pairing Hamiltonian 
that occurs when there is no translation or time-reversal 
symmetry restricting the states which can be paired. If 
<j) is small and we can use perturbation theory, the varia- 
tional energy difference between the superconducting and 
normal state, (H)s — (H)n, that follows from Eqs. (||) 
and (Q) is, keeping the lowest two orders in <fi: 
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This describes a continuous phase transition at the field 
H at which the quantity in parenthesis vanishes. When 
the electron gas wavefunctions are used, this yields the 
semiclassical result for H C 2 plus some cutoff-dependent 
corrections |J. 

Non-perturbative effects are illustrated most simply for 
a two-dimensional electron gas at zero temperature. Bare 
energies are Landau levels, multiples of the cyclotron en- 
ergy fioj H . The Fermi energy is at the highest occupied 
Landau level, so the matrix e in Eq. (Q) is 
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Near the phase transition 4> can be taken to be small 
compared to hw H , but before perturbation theory can 
be applied degenerate levels must be diagonalized. The 
pairing splits each Landau level into two with energies 
e m ± 4>\A m ^2Nf-m\- For the filled levels below the Fermi 
energy the number of occupied states shifted up in energy 
equals the number shifted down in energy so there is no 
net contribution to the energy from this splitting. When 
the level right at the Fermi energy is partially filled, how- 
ever, more states will be shifted down in energy (Fig. [I]). 
This will contribute a term linear in <j) to the energy dif- 
ference between the superconducting and normal state, 
which then has the form: 

AE 2D (H) = -e<f> + a{H)cf> 2 + 0(<P 3 ) . (10) 

In the relevant weak field limit, E F ^S> TiuJ H , and as- 
suming the lowest Landau level component for the or- 
der parameter, a calculation for the electron gas yields 
e = a pq fiLo H and a(H) = po ln(H/H® 2 ), where po is the 
zero-field density of states at the Fermi energy and H® 2 is 
the perturbative result for the upper critical field. Here 
a ~ 0.67 v/Ny 4 where v is the fraction of the highest 
Landau level which is filled (or one minus this fraction if 
it is more than half-filled) and Nf is the total number of 
filled levels N f = E F /hcu H = (k F £ H ) 2 /2 ~ 10 3 -10 4 for 
typical materials. 

Minimizing Eq. (|l^) for the energy difference yields 
the result shown in Fig. ^: the order parameter develops 
a "tail" of persisting superconductivity. In the limit of 
small e, the magnitude of the order parameter in the tail 
is (f> — e/2a(H). Since e depends on the filling factor 
of the highest level, cf> oscillates with magnetic field at a 
frequency 5H/H C 2 ~ l/Nf about the upper dashed line 
in Fig. |. 



This form of residual superconductivity depends on the 
exact degeneracy of the Landau level poles of a non- 
interacting two-dimensional electron gas in a magnetic 
field. In real materials interactions will broaden these 
levels. When the broadening is larger than ~ e/a(H), 
there will no longer be a term linear in 4> in the energy 
difference and this effect will be suppressed. 

It is believed that the effect of impurities in two di- 
mensions at zero temperature is to cause extended states 
to appear only at discrete energies, as evidenced by 
the quantum Hall effect (?]||. To the extent to which 
the most relevant interactions are those between the ex- 
tended states, pairing will split these levels as in Fig. 
and cause the two-dimensional residual superconduc- 
tivity to occur in superconductors at temperatures low 
enough and materials clean enough that the quantum 
Hall effect could be resolved in the normal state. Such 
conditions have not yet been achieved in thin metallic 
films, but may eventually be. 

Although the two-dimensional form of residual super- 
conductivity depends on large degeneracies in the spec- 
trum, a weaker form occurs more generally, when the 
normal-state spectrum is flat. To see where this comes 
from, consider first the zero-field BCS theory. The dif- 
ference between the ground state energy of the super- 
conductor and of the normal state can be written, for a 
uniform gap A, 

^E BCS = ^ - 2 po J C dE ( VE 2 + A2 - E) . (11) 

The first term is from the expectation value of the in- 
teraction term in H and the second is the shift in the 
electronic energy levels. The normal state spectrum is 
assumed flat in the range E c around the Fermi energy. 
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FIG. 1. Top: Density of states of an ideal two-dimensional 
electron gas. At zero temperature the Landau level at the 
Fermi energy is partially filled. Bottom: The pairing splits 
each level into two levels separated by an energy oc <j>, the 
magnitude of the order parameter. Electrons in the partially 
filled level contribute a net 0(<f>) term to the energy. 



FIG. 2. Solid line is the standard perturbative result for 
the order parameter as a function of field near the transition: 
(f> oc (H C 2 — H) 1 / 2 . Two-dimensional residual superconductiv- 
ity at T — 0, Eq. (p"(i|), yields rapid oscillations about the top 
dashed line; the three-dimensional form, Eq. ([l2]), yields the 
middle dashed line. 
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Even for small A compared to E c , perturbation theory 
can not be used to evaluate the integral in the second 
term of Eq. (|ll|). The exact answer is § A 2 [ ln(2^ c /A) + 
\ + 0(A 2 /E 2 )]: the A 2 In A term is non-perturbative. 
In the BCS case this term exactly cancels the A 2 /Vo 
term, leaving the net condensation energy ^poA 2 . In 
a magnetic field, the balance between the 0(</> 2 ) terms 
determines the perturbative upper critical field, as in Eq. 
(Q), and a net 4> 2 \sx<f> term remains: 

AE m (H) = e<p 2 \n<j) + a{H) </? + O(0 3 ) . (12) 

This form for the energy difference can be derived for 
the three-dimensional electron gas in the formal limit, as 
before, <p <C fiuj H . The additional complication is that 
the bare Hamiltonian does not have the simple structure 
of Eq. (Q) ; the energies are offset by an amount S z which 
depends on the ^-momentum: e n (k z ) = hui H (n + 5) + 
k 2 /2m—E F = huj H (n—Nf+5 z ). Integrating with respect 
to S z those levels within Tiu H of the Fermi energy gives 
an integral of the form which occurs in Eq. (pd|). The 
result is Eq. (|l2|) with e « 0.22 po/N 1 / 2 and a(H) = 
po HH/H° c2 ). 

The form of Eq. ( |l2| ) is not sensitive to the assump- 
tion of a non-interacting electron gas in a magnetic field. 
If we suppose that 1) the spectrum in the normal state 
is flat, and 2) the radiply varying matrix elements that 
enter Eq. (Q) are random, we obtain the results shown 
in Fig. ^. The diagonal matrix e is taken from a uni- 
form distrubition —N/2 < £j < N/2, corresponding to 
the flat normal state spectrum. The N x N matrix A a ^ 
is taken to be a complex symmetric matrix normalized 
to TrAty4 = N with individual elements Gaussian dis- 
tributed. (The 4> — > results are not very sensitive to 
this choice of distribution, or to a more realistic approxi- 
mation where A a/ 3 falls off exponentially with \e a — £p\.) 
We calculate the shift in the ground state energy AE, di- 
vide by (j) 2 , and examine the limit as 4> — > 0. If the result 
were purely perturbative, corresponding to Eq. (H), we 
would see a constant limit. Instead, a logarithmic upturn 
is clearly visible, corresponding to Eq. (|l2|). 

Whether the non-perturbative effects will be visible 
in a given material depends on the number of states N 
that are mixed by the pairing. In the electron gas case, 
the symmetry of the vortex lattice generates a conserved 
momentum and N ~ U Nf. The effects of disorder in a 
superconductor, when time-reversal symmetry has been 
explicitly broken by the presence of a magnetic field, 
are not well understood. (The semiclassical Abrikosov- 
Gor'kov description is not sufficient to address the effects 
discussed here.) If the number of states being paired be- 
comes macroscopic then non-perturbative effects will be 
suppressed. If the effective number of states being paired 
remains small, the number of states inside a volume de- 
termined by the scattering length rather than the co- 
herence length, then non-perturbative effects will remain 



and be enhanced. 

A further interesting possibiltity is that the random 
matrix model described above is a useful model for the 
spectrum of a granular superconductor || in a field, 
where pairing takes place within a grain and is in the 
mesoscopic regime. In that case N is the number of states 
in a grain being paired, and can be as low as 10 2 . 

At fields far above the conventional H C 2, minimiz- 
ing Eq. (|l^) for the order parameter strength gives 
4> ~ exp(— a(H)/e). The order parameter becomes ex- 
ponentially suppressed, as shown in Fig. ^. When <$> 
is sufficiently small the superconductivity will become 
disfavored due to normal-state Pauli paramagnetism or 
pair-breaking effects, such as magnetic impurities, which 
are not included in the idealized Hamiltonian in Eq. (jl]) . 
It is an interesting possiblity, however, that more exotic 
order parameter configurations [ fl0| would allow the resid- 
ual state to persist even in the presence of such effects. 

At finite temperature, the quantity to be minimized is 
the Gibbs free energy. As in the BCS case, the energy dif- 
ference AE(T) = <p 2 /V - 2 J dE E 6p{E)t&nh(E/2T), 
where Sp is the change in the density of states, and the 
entropy is given by the Fermi gas expression. The factor 
of tanh(£/2T) cuts off the (f> 2 /E 2 rise of 5p{E) at low 
energies which generated the 4> 2 In <j> term. Initial numer- 
ical results indicate that non-perturbative effects remain 
in the 0(<fi 3 ) contribution to the free energy. This would 
have interesting implications for behavior near the phase 
transition. 

Usually an O(cf) 3 ) term in the free energy 

G = ate)) 2 + 7 3 + b(j) 4 + ... , (13) 

where t = T — T C (H), signals a first order phase transi- 
tion. Defining 
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FIG. 3. Change in ground state energy vs. order parameter 
strength <f) for the random matrix ensemble defined in the text. 
Horizontal dashed line is the perturbative N — » 00 result. 
Points are for N = 5, 10, 20, (o, □, A) and solid lines are fits 
to the form of Eq. (jl|): AE/</> 2 = a log <p - j3 + 7 <j>. Both 
a and 7 decrease as ~ 1/N. 10 4 matrices were used for each 
point. 
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and keeping terms up to 0(<j> ), we have that a second 
minimum appears in G(cj>) for temperatures below t* and 
becomes energetically favorable compared to the = 
solution when t < t\ = (8/9)t*. More precisely, the 
two non-zero extrema of G which occur at temperatures 
below t* are at 



-l±Vl -(*/** 
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For negative 7, + is the local minimum which becomes 
energetically favorable below t±, causing a first-order 
transition. For positive 7, the local minimum <f>_ < 0. 
0, however, is the magnitude of the order parameter and 
must be positive, so this solution is not allowed. What 
happens instead, if 7 is positive, is that there is a con- 
tinuous transition at t — as the 0+ solution becomes 
an allowed minimum. Near this transition, for \t\ -C t* , 
expanding the square root yields 



30* 



(16) 



The scale for this deviation from the familiar mean-field 
behavior oc is set by t*: for |f| 3> t* the square- 
root behavior is recovered. The same behavior occurs if 
the field is varied at fixed temperature; t is replaced by 
h = H - H c2 (T) and the result is oc \h\ for \h\ < h* . 
Although the models considered here give a negative 7, 
favoring the continuous transition, the possibility that 
additional physics enters at sufficiently small 0, and gen- 
erates a positive 7, can not be ruled out. 

The deviations from standard mean field behavior de- 
scribed here would be difficult to see in resistivity mea- 
surements, the most common measurements done near 
H C 2 on conventional materials. They would be more 
likely to be seen in sensitive heat capacity, magnetiza- 
tion, or tunneling measurements. The form Eq. ( |l3|) for 
the free energy gives G oc (T — T c ) 3 near the continuous 
transition (for t <C t*), instead of the Ginzburg-Landau 
behavior G oc (T — T c ) 2 , and leads to a kink in the heat 
capacity rather than an abrupt discontinuity. The mag- 
netization M oc 2 would vary as (T c — T) 2 rather than 
(T c — T). In the perturbative mean field theory, the spec- 
trum near H C 2 changes by O(0 2 ). Non-perturbative ef- 
fects arise from states close to the Fermi energy being 
strongly affected by the pairing even as — > 0. This 
leads to a zero-bias anomaly in the spectrum: a peak 
near the Fermi energy, a (fx 2 / E 2 rise, at the center of 
which is is a sharp spike giving p(E F ) = 0. 

Non-perturbative effects appear at scales much larger 
than the deviations from mean field behavior expected 
from phase fluctuations. The phase fluctuations which 
cause vortex lattice melting are believed to become im- 
portant on the scale of critical phenomena, as esti- 
mated by the Ginzburg criterion [0. This scale is 



(Tq - T c )/T c < 1(T 16 in zero magnetic field and < 1CT 7 
at H = 10 T for typical parameters of bulk conventional 
superconductors, although the scale may be much larger 
in high-T c materials Oj. The scale for the deviations 
described here is 1/ y/Nf which can be 1%— 10% for con- 
ventional supercondcutors. 
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